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Abstract 

We calculate the branching ratios and CP-violating asymmetries for B° — » K°K , K K*°, 
K + K*~, K~K* + , and B + — ► K + TC* and K K* + decays by employing the low energy ef- 
fective Hamiltonian and the perturbative QCD (pQCD) factorization approach. The the- 
oretical predictions for the branching ratios are Br(B°/B — ► K^K*^) ~ 7.4 x 10~ 8 , 
Br{B°£B° -► K K*°(K°K* )) « 19.6 x 10~ 7 , Br(5+ -> K+iT°) « 3 x 10" 7 and Br(5+ -► 
) « 18.3 x 10~ 7 , which are consistent with currently available experimental upper lim- 
its. We also predict large CP-violating asymmetries in these decays: Agp^K^K ) ~ —20%, 



A^piK^K ) « -49%, which can be tested by the forthcoming B meson experiments. 
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I. INTRODUCTION 



The study of exclusive non-leptonic weak decays of B mesons provides not only good 
opportunities for testing the Standard Model (SM) but also powerful means for probing 
different new physics scenarios beyond the SM . The mechanism of two body B decay is still 
not quite clear, although many scientists devote to this field. Starting from factorization 
hypothesis [l|, many approaches have been built to explain the existing data and some 
progresses have been made. For example the generalized factorization (GF)|2J, QCD 
factorization (QCDF) approach [1,0], the perturbative QCD (pQCD) approach p, @, 0, [1] 
and the soft-collinear effective theory (SCET) [§]. The pQCD approach is based on Kt 
factorization theorem [To|] while others are mostly based on collinear factorization fl~l| . 

In our opinion, the pQCD factorization approach has three special features: (a) Su- 
dakov factor and threshold resummation [jjj are included to regulate the end-point singu- 
larities, so the arbitrary cutoff [l3[ is no longer necessary; (b) the form factors for B — > M 
transition can be calculated perturbatively, although some controversies still exist about 
this point; and (c) the annihilation diagrams are calculable and play an important role 
in producing CP violation T^. Up to now, many B meson decay channels have been 



studied by employing the pQCD approach, and it has become one of the most popular 
methods to calculate the hadronic matrix elements. 

In this paper, we will study the branching ratios and CP asymmetries of B — > K K* 
decays in the pQCD factorization approach. Theoretically, in the B — > KK* decay 
modes, the B meson is heavy and sitting at rest. It decays into two light mesons with 
large momenta, so these two energetic final state mesons may have no enough time to 
get involved in soft final state interaction (FSI). In this case, the short distance hard 
process dominates the decay amplitude and the non-perturbative FSI effects may not be 
important, this makes the pQCD approach applicable. At the same time, the B — ► K K* 
decays have been studied before in the GF approach [2;] and the QCDF approach [!, 0]. 
The similar decays such as B — > KK and K*K* decays have been investigated in the 
pQCD approach recently [l5|, T6[. On the experimental side, the first measurement of 



B» -> (K°K + TK* ) decay has been reported very recently by BaBar collaboration 



17j | in units of 10 6 (upper limits at 90% C.L.' 

Br(B° - K°1T + K°K*°) = 0.2^+^ (< 1.9). (1) 



-*Q 



For B + — * K + K decay, only the experimental upper limit is available now [18|, Il9 



T*0 



Br(B + -> K + K™) < 5.3 x 10~ 6 . (2) 



This paper is organized as follows. In Sec. II, we give the theoretical framework of the 
pQCD factorization approach. Next, we calculate the relevant Feynman diagrams and 
present the various decay amplitudes for B — > KK* decays. In Sec. IV, we show the 
numerical results of the CP-averaged branching ratios and CP asymmetries and compare 
them with currently available experimental measurements or the theoretical predictions 
in QCDF approach. The summary and some discussions are included in the final section. 
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FIG. 1: Factorization for B — > KA"* Decays 



II. THEORETICAL FRAMEWORK 

The three scales pQCD factorization approach 0, 0] has been developed and applied 
in the non-leptonic B meson decays for some time. In this approach, the decay ampli- 
tude is factorized into the convolution of the mesons' light-cone wave functions, the hard 
scattering kernel and the Wilson coefficients, as illustrated schematically by Fig. 1, which 
stands for the soft, hard and harder dynamics characterized by three different energy 
scales (t 



0(-\/ AM B ),m b , M w ) respectively. Then the decay amplitude A(B 
is conceptually written as the convolution 



M X M 2 



A(B -> M X M 2 ) ~ Jd%d 4 k 2 d% Tr [C(t)®B(ki)<&M 1 (k2)®M 2 (k3)H(ki, h, k 3 ,t)] , (3) 

where fcj's are momenta of light quarks included in each mesons, and the term "Tr" 
denotes the trace over Dirac and color indices. C(t) is the Wilson coefficient which 
results from the radiative corrections at short distance. In the above convolution, C(t) 
includes the harder dynamics at scale larger than M B and describes the evolution of 
local 4-Fermi operators from mw (the W boson mass) down to t ~ scale, 
where A = M B — m b . The function H(ki, k 2 , k 3 , t) describes the four quark operator and 
the spectator quark connected by a hard gluon whose q 2 is of the order of AM B , and 
includes the 0(^/AM B ) hard dynamics. Therefore, this hard part H can be evaluated 
as an enpansion in power of as(t) and A/t, and depends on the processes considered. 
The function $^ {M = B, Mi, M 2 ) is the wave function which describes hadronization 
of the quark and anti-quark into the meson M, and independent of the specific processes. 
Using the wave functions determined from other well measured processes, one can make 
quantitative predictions here. 

Since the b quark is rather heavy we consider the B meson at rest for simplicity. It is 
convenient to use light-cone coordinate {p + ,p~,Px) to describe the meson's momenta, 



V 



± 



— (p°±p 3 ), and PT = <y,p 2 ) 



(4) 



Using the light-cone coordinates the B meson and the two final state meson momenta can 
be written as 



Pi = ^§(M,Ot) 



M, 



V2 



z?/ l,r|,,0 T ) 



M 



B 



V2 



(0,l-ri,O T ) 



(5) 
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respectively, where tk* = m^/n^B] and the terms proportional to m 2 K /m B have been 
neglected. 

For the B — > JCK* decays considered here, only the K* meson's longitudinal part 
contributes to the decays, its polarization vector is €l = (1, — r^*, T ). Putting the 

light (anti-) quark momenta in B, K* and K mesons as k\, k 2 , and k 3 , respectively, we 
can choose 



h = (xiP+, 0, k 1T ), k 2 = (x 2 P+, 0, k 2T ), k 3 = (0, x 3 P 3 ~, k 3T ) 
Then the integration over k± , k 2 , and /cj|" in Eq.Q will lead to 

A(B — >• KK*) - 



(6) 



dx\dx2dx 3 b\db]b2db2b 3 db 3 

Tr [C(t)^ B ( Xl , 6i)$ fc *(x 2 , 6 2 )$ k (x 3 , h)H(xi, b h t)S t (x t ) e" 5(i) ] ,(7) 

where bi is the conjugate space coordinate of k^, and t is the largest energy scale in 
function H(xi,bi,t). The large logarithms hi(mw/t) coming from QCD radiative cor- 
rections to four quark operators are included in the Wilson coefficients C(t). The large 
double logarithms (ln 2 Xj) on the longitudinal direction are summed by the threshold 
resummation l3] , and they lead to S t (xi) which smears the end-point singularities on 
Xj. The last term, e~ s<yt \ is the Sudakov form factor resulting from overlap of soft and 
collinear divergences, which suppresses the soft dynamics effectively[20]. Thus it makes 
the perturbative calculation of the hard part H applicable at intermediate scale, i.e., Mb 
scale. 

The weak effective Hamiltonian H e ff for B — > KK* decays can be written as 21 



Keff ~ — 



G f 

71 
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v ub v: d (Ci(/i)oy(A*) + c 2 (fi)o u M) - v tb v; d £ c^) o^) 



i=3 



(8) 



where Cj(/i) are Wilson coefficients evaluated at the renormalization scale fi and Oi are 
the four-fermion operators for b — > d transition: 
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o 7 
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o 8 = 


o 9 


= \d a YLb a 


Oio = 



da^LUa ■ UfjJpLbp , 

daYLbp ■ Eg' q'pitiRq'a > 
\d_ a YLb p ■ 2g' e q >q'pinRq' a 
\d a YLb p ■ J2 q > Zq'q'piM'a 



(9) 



where a and /3 are the SU(3) color indices; L and R are the left- and right-handed 
projection operators with L = (1 — 75), -R = (1 + 75). The sum over runs over the quark 
fields that are active at the scale \i = 0(m b ), i.e., (q'e{u, d, s, c, b}). For the decays with 
b — > s transition, simply make a replacement of d by s in Eqs. (JSJ) and ([9]). 

The pQCD approach works well for the leading twist approximation and leading double 
logarithm summation. For the Wilson coefficients Ci(fi) (i = 1,...,10), we will also 
use the leading order (LO) expressions, although the next-to-leading order calculations 



already exist in the literature [21]. This is the consistent way to cancel the explicit /x 



dependence in the theoretical formulae. For the renormalization group evolution of the 



4 



Wilson coefficients from higher scale to lower scale, we use the leading logarithmic running 



equations as given in Appendix C and D of Ref. [22 



In the resummation procedures, the B meson is treated as a heavy-light system. In 
general, the B meson light-cone matrix element can be decomposed as [23| 



1 d 4 z 



i 



iki z 



(0\b a (0)dp(z)\B{p B )) 



V2N C 



>B 



's(ki 



ifi + — 



■0s(ki 



(10) 



/3a 



where n + = (1, 0, T ), and n_ = (0, 1, T ) are the unit vectors pointing to the plus and 
minus directions, respectively. From the above equation, one can see that there are two 
Lorentz structures in the B meson distribution amplitudes. They obey to the following 
normalization conditions 



(2tt) 4 



>B 



Mi) 



B 



2y/2N c 



d 4 h 

WT 4 



0. 



11) 



In general, one should consider these two Lorentz structures in calculations of B meson 
decays. However, it can be argued that the contribution of 4>b is numerically small 24j . 



thus its contribution can be numerically neglected safely. Using this approximation, we 
can reduce one input parameter in our calculation. Therefore, we only consider the 
contribution of Lorentz structure 



1 



B 



V2Nc 



h + "ifi)750i?(ki 



(12) 



The K and K* mesons are treated as a light-light system. Based on the SU(3) flavor 
symmetry, we assume that the wave functions of K and K* mesons are the same in 
structure as the wave functions of 7r and p, respectively, then the K meson wave function 



is defined as 25, 26 



:75 



K 



,K±P , 



K 



vn)4> T K (x)} 



(13) 



where P and x are the momentum and the momentum fraction of K, respectively. The 
parameter ( is either +1 or —1 depending on the assignment of the momentum fraction x. 
While in B — > KK* decays, K* meson is longitudinally polarized, only the longitudinal 
component of the wave function should be considered [H, 23], 



<T) L 

K* 



V2Nc 



{{ [$k*4>k*{x) + m K *(f) K *(x)] +m K *(j) S K - t (x)} 



(14) 



The second term in above equation is the leading twist wave function (twist-2), while the 
first and third terms are sub-leading twist (twist-3) wave functions. The transverse part 
of can be found for example in Ref. [161 ] . 

The explicit expressions of the distribution functions B (kx), ^(x), 4>x(x), ^(x), 
4>k*{x), 4>k*{ x ) an d 4>k*{x) will be given in next section. The initial conditions of leading 
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Mi = K°(K*°) 
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2^ 



(c) 
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FIG. 2: Typical Feynman diagrams contributing to B° — > if ^* ^* .^ °) decays. The diagram 
(a) and (b) contribute to the form factor Aq^ k * or i^-j - ^ for Mi = if* or if , respectively. 
Other four Feynman diagrams obtained by connecting the gluon lines to the d quark line inside 
the B° meson for (e) and (f), and to the lower s or d quark line for (g) and (h) are omitted. 

twist distribution functions 4>i(x), i = B, K*, K, are of non-perturbative origin, satisfying 
the normalization condition 



(j>i(x,b = 0)dx = T^fi , 
where /j is the decay constant of the corresponding meson. 



(15) 



III. PERTURB ATIVE CALCULATIONS 



For the considered decay modes, the Feynman diagrams are shown in Figs. [2H We 
firstly analyze the corresponding decay modes topologically: (i) the eight diagrams can 
be categorized into emission and annihilation diagrams; (ii) each category contains four 
diagrams: two factorizable and two nonfactorizable. In Fig. 2, for example, Figs. 2(a-d) 
are emission diagrams, while Figs. 2(e-h) are annihilation ones topologically; and Figs. 
2(a,b,g,h) are factorizable and Figs. 2(c-f) are nonfactorizable diagrams. 

For B° — > K°K (K*°K ) decays, only the operators O^iq contribute via penguin 
topology with light quark q = s (diagrams a,b,c,d ) and via the annihilation topology 
with the light quark q = d (diagram 2(f) and 2(h) ) or s ( diagram 2(e) and 2(g) ). It is 
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FIG. 3: Feynman diagrams for B° —> K + K* (K* + K ) decays. 

a pure penguin mode with only one kind of CKM elements, and consequently, there is no 
CP violation for these decays. 

For the B°(B°) — > K + K*~ (K* + K~) decays ( see Fig. [3D, the current-current operators 

0^1 contribute via the annihilation topology ( Figs. 3(c,d,g,h)), while the operators O3„io 
contribute via the annihilation topology with the light quark q — s ( Figs. 3(a,b,e,f)) or 
q = u ( Figs. 3(c,d,g,h)). 

For the B + — > K + K*°(K* + K°) decays ( see Fig. H}, the current- current operators O^l 
contribute via the annihilation topology ( Figs. 4(e-h)), while the penguin operators O 3 _i 
contribute via the penguin topology with the light quark q = s (Figs. 4(a-d)) or via the 
annihilation topology with q = u ( Figs. 4(e-h)). 

In the analytic calculations, the operators with (V — A)(V — A) structure work directly, 
while the operators with (V — A)(V + A) structure will work in two different ways: 

• In some decay channels, some of these operators contribute directly to the decay 
amplitude in a factorizable way. 

• In some other cases, we need to do Fierz transformation for these operators to 
get right flavor and color structure for factorization to work. In this case, we get 
(S + P)(S- P) operators from (V - A)(V + A) ones. 
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B+ 




FIG. 4: Feynman diagrams for B + — > K + K*°(K* + K ) decays. 
A. 5° -> K K*°(K* K°) decay 

For the sake of the reader, we take the S° — > A 1 ^* ^* !^) decay channel as an 
example to show the ways to derive the decay amplitude from individual diagram. As 
shown explicitly in Fig. 2(a), the meson Mi which picks up the spectator quark can be 

A or A* , the emitted meson M 2 should be A" or A at the same time. The B meson 
therefore can decay into the final state / = A°A and / = A*°A simultaneously. 
The B meson, on the other hand, also decay into the same final state / = A A* and 
/ = A*°A simultaneously. 

Now we consider the usual factorizable diagram 2(a) and 2(b) for the case of M\ = A* . 
The (V — A)(V — A) operators 0^,a and Og^o contribute through diagram 2(a) and 2(b), 
the sum of their contributions is given as 

F e K* = ^y/^GF^CpfK^B / dxidx 3 / bidbib s db s (f)B(xi,bi) 

Jo Jo 

■{[(l + x 3 )4> K *{x 3 ,b 3 ) + (l-2x 3 )r K *((f) s K *(x 3 ,b 3 ) + (t) t K ,{x 3 ,b 3 ))] 
■a s (tl)h e (xx, x 3 , bi, b 3 ) exp[-S a {tl)] 

+2r K *4> s K *(x 3 ,b 3 )a s (t 2 e )h e (x 3 ,x 1 ,b 3 ,b 1 )exp[-S a (t 2 e )}} , (16) 
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where CV = 4/3isa color factor. The functions h\, the scales t\ and the Sudakov factors 
S a (tl) and S a (tl) will be given explicitly in the Appendix. In Eq. f[T6"j) . we do not include 
the Wilson coefficients of the corresponding operators, which are process dependent. They 
will be shown later in this section for different decay channels. 

The form factor of B to K* transition, A B ^ K * (0), can also be extracted from F e x* in 
Eq. (USD, that is 

The operators 0$-$ have a structure of (V — A)(V + A). Some of these operators 
contribute to the decay amplitude in a factorizable way. Since only the axial- vector part 
of (V + A) current contribute to the pseudo-sealer meson production, 

(K*\V - A\B){K\V + A|0> = -(K*\V - A\B){K\V - A|0). (18) 

The contribution of these operators is opposite in sign with F e x* in Eq. (|T6l) : 

F e p K \ = -F eK *. (19) 

In some other cases, one needs to do Fierz transformation for these operators first and 
then get right color structure for factorization to work. In this case, one gets (S—P)(S+P) 
operators from (V — A)(V + A) ones. For these (S — P)(S + P) operators, Figs. 2(a) 
and 2(b) gives 

rl poo 

FeK* = 8\/2G fttC f fx r K m % / dx%dx 3 / b±dbib 3 db 3 </>b(xi, h) 

Jo Jo 

■ { [<pK*( x 3,h) + r K *{{x 3 + 2)(p s K ,(x 3 ,b 3 ) - x 3 $ K *{x 3 , 6 3 ))] 
■a s {t\)h e {x u x 3 , bi, b 3 ) exp[-S a (tl)) 

+ (xi(j)K* (%3, 63) + 2r K ,(j) s K , (x 3 , b 3 )) a s (t 2 e )h e (x 3 , x u b 3 , b{) exp[-S a (l%)] } (20) 

For the non-factorizable diagram 2(c) and 2(d), all three meson wave functions are 
involved. The integration of 63 can be performed using 5 function 5(b 3 — b±), leaving 
only integration of b\ and b 2 . M e x* denotes the contribution from the operators of type 
( V — A) ( V — A) , and is the contribution from the operators of type ( V — A) ( V + A) : 

16 f°° 
M eK * = —j=G pirC Vmg / dx\dx 2 dx 3 / b\db^b 2 db 2 0b(xi, &i)0^-(x 2 , b 2 ) 

v3 io Jo 

■ {- [-x 2 (j)K*{x 3 M) + r K *x 3 {(j) s K *{x 3 M) - 4>K*{x 3 ,bi))] 
■a s (t f )h 1 f (x 1 ,x 2 , x 3 , 61, b 2 ) exp[-5 c (t})] 

+ [(^2 -x 3 - l)(f) K *(x 3 , b x ) +r K *x 3 ((j) s K ,{x 3 M) +0^*(^3,&i))] 
■a s (t f )h 2 f (x 1 ,x 2 ,x 3 ,b 1 ,b 2 )exp[-S c (t 2 f )]} , (21) 
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16 f 1 f°° 

M^x* = —j=G p^C ftti% \ dxidx 2 dx 3 / bidbib 2 db 2 (f>B(xi,bi)r K 
v3 Jo Jo 

■ { [(xi - x 2 ) (<pK( x 2,b 2 ) ~ <pK(x 2 ,b 2 )) 4>k*{xzM) +r K * (xi (<pK( x 2,b 2 ) 
-<p T K (x 2 ,b 2 )) (^.(x 3 ,6i) - 0x-.(x 3 ,6i)) -X2 (<££(z2,&2) -^(^2,62)) 
• (^.(3:3,61) - ^.(x 3 ,6i)) - x 3 (0£(:r 2 ,6 2 ) + 0^2, ^2)) (0^*(x 3 ,&i) 
+0^.(^3, &i)))] a s (t/)/i}(xi, x 2 , x 3 , 61, 6 2 ) exp[-5' c (^)] 
- [(xi + x 2 - 1) (0^(x 2 , 62) - 0x(a: 2 , ^2)) 0x* (x 3 , h) 
+r K * (xi (<j)K(x 2 ,b 2 ) - <p T K (x 2 ,b 2 )) {(/>K*(x3,h) - 4>k*(x3, b i)) 
-(l-x 2 ) ((/)K(x2,b 2 ) - 4> T K{x 2 ,b 2 )) (<f> a K .(x 3 ,bi) - 0' x *(x 3 ,&i)) 

-x 3 {<P p K (x 2 ,b 2 ) + 4> T K (x 2 M)) {<t>K*{xzM) + 0^(^3,61)))] 

Q! i (t/)^(a:i,X2,a:3,6i,62)exp[-S , c (^)]} . (22) 

For the non-factorizable annihilation diagram 2(e), we have three kinds of contribu- 
tions: M aK * for (V - A)(V - A) operators, for (V-A)(V + A) operators and 
for (S -P)(S + P) operators. 

16 f 1 f°° 

M a K* = —^Gf^Cftti^ / dx\dx 2 dx 3 I bidbib 2 db 2 (f)B(xi,bi) 



V3 jo jo 

■ {- [x 2 (j)K*(x3,b 2 )(l)^(x 2 ,b 2 ) +r K *r K ((p P : (x 2 ,b 2 ) ((x 2 + x 3 + 2.) 

■<Pk* (x 3 , 62) + (X2 - x 3 ) (fix* (x 3 , 6 2 )) + <Pk( x 2, b 2 ) (-x 3 ((j) s K * (x 3 , b 2 ) 
-4>K* (*3, 62)) - 20^* (X 3 , 62) + ^2 (Xa, b 2 ) + 0*^ (X 3 , 62)) ) )] 

■a s (t 3 f )h 3 f (x 1 ,x 2 , x 3 , 61, 6 2 ) exp[-5' c (^)] 

+ [^30^*(^3,&2)0^(^2,&2) -r K *r K (-x 2 (0^(x 2 ,& 2 ) - 0^(x 2 , 6 2 )) 
•(0^-0^0 -a* (0^(x 2 ,6 2 ) +0^(x 2 ,6 2 )) (0^(x 3 ,fc 2 ) + 0^*(x 3 ,6 2 )))] 
•a s (t))/i)(a;i,x 2 ,X3,6i,6 2 )exp[-S' c (t))] } , (23) 

where tk = /ms with tUq = m 2 K /{m s + rrid). 

16 

^oif* = —j=GF^CFm A B I dx\dx 2 dx 3 I b 1 db i b 2 db 2 (j)B(x 1 ,bi) 
v3 Jo Jo 

■ { [-r K * (x 3 - 2) ((f> a K , (x 3 , 62) + (ar 3 , 62)) 

(x 2 - 2) 0^»(a; 3 ,6 2 ) (0^(x 2 ,6 2 ) + 0^(x 2 ,6 2 ))] 
•a a (l? f )h 3 f (xi, x 2 , x 3 , 61, & 2 ) expf-^c^J)] 
+ [~a; 2 r^0^*(x 3 ,6 2 ) (0^(x 2 ,6 2 ) + 0^(x 2 ,6 2 )) 
+x 3 r ft: »0^(a; 2 ,& 2 ) (0^*(x 3 ,6 2 ) + 0* K »(x 3 , 62))] 

•a s (t))/i)(a; 1 ,X2,a; 3 ,6i,6 2 )exp[-S' c (t))]} . (24) 
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16 f 1 f°° 

M^x* = —=G F irC F m 4 B / dx\dx 2 dx 3 / b\db]}o 2 db 2 <Ab(£i, b\) 
v3 Jo Jo 

■ { [x 3 (f)K*{x 3 ,b 2 )(f)^ : (x 2 , b 2 ) + r K *r K (((x 2 + x 3 + 2) 0^,(0:3,62) 

- ( x 2 -X 3 ) 0^.(2:3,62)) 0X^2,62) + (X3 (0^,(X 3 ,6 2 ) + 0X*( a; 3 ! &2)) 

+0:2 (0^(x 3 ,6 2 ) -cf) s K *(x 3 ,b 2 )) - 2(f) t K *(x 3 ,b 2 )) (f)J c (x 2 ,b 2 ))] 

■ a s (t 3 f )h 3 f (x 1 ,x 2 , x 3 , 61, 6 2 ) expf-S^tJ)] 

+ [-x 2 (j>K*{xz, b 2 )(f)x(x 2 , b 2 ) +r K *r K (-x 2 ((f)^(x 2 ,b 2 ) + <fr K (x 2 ,b 2 )) 

• {4>K*( x ^b 2 ) + (/> K *( x 3,b 2 )) - x 3 (0^(x 2 ,6 2 ) - 0^(^2,62)) 

• (0^*(^3, b 2 ) -0^(x 3 ,6 2 )))] a s (t))/i)(a; 1 ,X2,X3,6i,62)exp[- 1 S c (t))]} .(25) 

The factorizable annihilation diagram 2(g) involves only K* and wave functions. The 
decay amplitude F a K*, F^k* and F^k* represent the contributions from (V — A)(V — A) 
operators, (V — A)(V + A) operators and (S — P)(S + P) operators, respectively. 

pi poo 

FaK* = -4:V2TcG F C F f B m 4 B / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

■ { [x 3 ((>k* (x 3 , 63)0^(^2, b 2 ) + 2r K *r K (p^(x 2 , b 2 ) ((1 + x 3 )(f) s K * (x 3 , b 3 ) 
-(1 - x 3 )<l> t K ,(x 3 , b 2 ))] a s (t 3 e )h a (x 2 , x 3 , b 2 , b 3 ) exp[-S d (l%)] 

- [x 2 (pK*(x 3 ,b 3 )(j)^(x 2 , b 2 ) +2r K *r K (f) s K ,(x 3 ,b 3 ) ((1 + x 2 )<f)^(x 2 ,b 2 ) 

-(1 - x 2 )0^(a;2, b 2 ))] a s (t 4 e )h a (x 3 , x 2 , 63, b 2 ) exp[-S d (t 4 e )] } , (26) 

F aK* = - F aK*, (27) 

n 1 /»oo 

F aK* = -8V2G F 7rC F m 4 B f B / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

■ { [2r K (f) K *(x 3 ,b 3 )(f)K(x 2 ,b 2 ) + x 3 r K * (<f) s K *(x 3 ,b 3 ) - ^(x^b^) 0^(x 2 ,6 2 )] 
•a s (t 3 e )h a (x 2 , x 3 , b 2 , b 3 ) exp[-5' d (^)] 

+ [2r K *(p s K ,(x 3 ,b 3 )(p^(x 2 ,b 2 ) + x 2 r K (0^(x 2 ,6 2 ) - <P K ( x 2,b 2 )) <pK*(x 3 ,b 3 )] 
■a s (t 4 e )h a (x 3 ,x 2 ,b 3 ,b 2 )eM-S d (ti)}} ■ (28) 

In the above equations, we have assumed that X\ « x 2 ,x 3 . Since the light quark mo- 
mentum fraction X\ in B meson is peaked at the small X\ region, while quark momentum 
fraction x 2 of K is peaked around 0.5, this is not a bad approximation. The numerical 
results also show that this approximation makes very little difference in the final result. 
After using this approximation, all the diagrams are functions of k\ = X\m B j\[2 of B 
meson only, independent of the variable of k±. 

For the Feynman diagram 2(f) and 2(h), the corresponding decay amplitude is the same 

in structure as those for 2(e) and 2(g). We get the decay amplitude easily by making two 

replacements of x 2 — > 1 — x 2 and x 3 — > 1 — x 3 in the relevant distribution amplitudes. 

*o 

For the case of Mi = K and M 2 = K , by following the same procedure, one can 
find all decay amplitudes: F eK , F^ and F e % M eK , M% M aK , M* and M% F aK , F^ K 
and F^K. The explicit expressions of these decay amplitudes will be given in Appendix 
A. 



11 



B. Total Decay amplitudes 



Based the isospin symmetry and the analytical results obtained in last subsection, 
one can derive out all the decay amplitudes for B° — > K + K*~ (K* + K~) and B + — > 
K+K*°(K* + K°) decays. 

Combining all contributions, the total decay amplitude for all considered decay modes 
can be written as: 

M(B° K«K-°) = -(, \f, k + M, K (ft - f ) + <i (ft - f ) 

+M, K (ft + ft - £ - %i) + M« (ft - 2i) + M« (ft - f ) + M„ K . (ft - %>) 

+ ^(ft + f-ft-T + T + T-T-^) + ^'(f + ft-f-f)} < 29 > 



M(B° K-Tf) = -£, (^ + ft - ^ - ^) + F5. (^ + ft - 2l - ^ ) 

+M. K . (ft - 5S) + M«. (ft - f ) + (ft + ft " f " 

+ M&. (ft - f ) + Mft. (ft - f) + *ft. + ft - ^ ~ °i) 
+ ^.(ft + |-ft-^ + ^ + f-f-£p) + M. K (ft-^) 

+f„ (5ft + 5 C « - C = - T + f + T - 5 C « - 5 C '») + (ft - t) } < 30 > 



+F aK 
+F aK * I < 



M a ^C 2 + F aX ( Ci + ^ 



n , C4 r , C 6 „ C 8 C10 
C 3 + T - C 5 - T - C7 7 - T + C7 9 + — 



6 |M aK (C 4 + C10) + M^.. (c 6 - 



3 



C 5 



3 2 



6 



1 



:C 9 



C10 
~6~ 



+ MaK* C 4 



C10 



(31) 



-» K*+R- 



(7 ^6 | C*7 | Cs Cg 



C10 
~6~ 



Wftft (c 3 + ^-C 5 -^-C 7 -^ + C 9 + ^) ft .w,„> ;r r M,„i 



+M oif C 4 



( -f) ^(c 6 ( -t) K«-r h) 



(32) 
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M{B A 



K+K*°) = ( M aK d 



-t.t\F eK [— + C A -- - 



10 



+F. 



° 5 -i- r j-° 7 



Ok — 



+ (C, c„ i + if + rvi + /; K ( ^ + c A + ^ + c 10 



(33) 



M (B + -> ) = & ( M^.d + F aA - 



p 2 

A'* 



C5 
3 

C 5 



C 6 



C 7 



C7 Cg 
~6 2 



F P 2 



c. 



c 7 



3 



C 4 



C 3 



C9 _ C10 
2 



+ M aK , (C 3 + C 9 ) + M^, (C s + C r ) + FaA 



where £ u = V* h V u d, £t — ^Vtd- The exact expressions of individual transition amplitudes 
not given explicitly in this section, such as F aK and M aK , etc., are collected in the 
Appendix A. 

The decay amplitudes for those charge-conjugated decay channels can be obtained from 
the results as given in Eqs. (l29|) - (J34J) by simple replacements of £ u — > £* and £t — > Q. 

Analogous to Eq. dill), the form factor F^ K (q 2 = 0) can also be extracted from F e K 
via the following relation 



V2F ( 



e.K 



G F f K *m 2 B 



(35) 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



A. Input parameters and wave functions 

Before we calculate the branching ratios and CP violating asymmetries for the B de- 
cays under study, we firstly present the input parameters to be used in the numerical 
calculations. 



(/=4) 



A— 

MS 



0.25GeV, f B = 0.19GeV, 



m, 



K 



1.7GeV, 



f K * = 0.217GeV, f K « = f K = 0.16GeV, m K = 0.497GeV, 
m K * = 0.89GeV, M B = 5.2792GeV, M w = 80.41GeV. (36) 

The central values of the CKM matrix elements to be used in numerical calculations are 



\V ud \ 
\Vtb\ 



0.9745, 
0.9990, 



\V ub \ 
\Vtd\ 



0.0036, 
0.0075. 



For the B meson wave function, we adopt the model 15j, [22j, |24 

4>B(x,b) = Nbx 2 (1 — x) 2 exp 



M 2 X 2 l f ,.r 



2u 2 



(37) 



(3* 
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where the shape parameter ojj, = 0.4 ± 0.04 GeV has been constrained in other decay 
modes. The normalization constant Ng = 91.745 is related to fs = 0.19GeV and uj^ = 0.4. 

The K* meson distribution amplitude up to twist-3 are given by {27] with QCD sum 
rules. 



1 + 0.57(1 - 2x) + 0.07C 2 3/2 (1 - 2x) 



JK * {0.3(1 - 2x) (3(1 - 2x) 2 + 10(1 - 2x) - l) 



(39) 



J K* 



X 



2VQ 

+ 1.68C 4 1/2 (1 - 2x) + 0.06(1 - 2xf (5(1 - 2xf - 3) 
+0.36 [1 - 2(1 - 2x) - 2(1 - 2x) ln(l - x)}} , 

fx* 



2V6 



{3(1 - 2x) [1 + 0.2(1 - 2x) + 0.6(10x 2 -10x+l 



-0.12x(l -x) + 0.36 [1 - 6x - 2 ln(l - x)}} , 
where the Gegenbauer polynomials are defined by 

Cl'\t) = 5 (5t 2 - 1) , C\ /2 {t) = i (35t 4 - 30t 2 + 3) . 
Z o 



(40) 



(41) 



(42) 



For K meson, we use 0^ of twist-2 wave function and (f)^ and 4>k of the twist-3 



wavef unctions from [26|, |27 

3 



x 



(f) K (x) 



4> T K {x) 



y/E 
Sk 



2V6 



f K x(l -x)[l + 0.51(1 - 2x) + 0.3 (5(1 - 2a;) 2 - l)] , 
[1 + 0.12(3(1- 2a;) 2 - 1) 



-0.12 (3 - 30(1 - 2x) 2 + 35(1 - 2x) 4 ) /8] 



^=(1 - 2x ) i 1 + °- 35 ( 10a;2 - 10x + !)] • 



(43) 

(44) 
(45) 



Based on the definition of the form factor Aq and Fq{* as given in Eqs. (fl7|) and 
(|35|) . we find the numerical values of the corresponding form factors at zero momentum 
transfer. 



-4, 



B^K* 



(q 2 = 0) 



0.46+0-07, 



-0.06 



Mb), 



(46) 



where the errors are induced by the change of cj& for Up — 0.40 ± 0.04 GeV. These results 
are close to the light-cone QCD sum rule predictions [28] 



Ar K \q 2 = 0) 

? B-+K/2 



0.374 + 0.034, 
0.331 + 0.041 



(47) 
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B. Branching ratios 



In order to calculate the branching ratios and CP asymmetries in a more clear way, we 
rewrite the decay amplitudes as given in Eqs. (l29l) - (l34"l) in a new form 

(48) 



M = V* b V ud T - V*V td P = V* b V ud T [I + ze^} , 



where the term "T" and "P" denote the "Tree" and "Penguin" part of a given decay 
amplitude M., which is proportional to £ u = V* b V ud or £ t = V t * b Vt d , respectively While 
the ratio 



(49) 



v; b v td 




p 


v* b v ud 




T 



is proportional to the ratio of penguin (P) to tree (T) contributions, the CKM angle 
is the weak phase, and S is the relative strong phase between the tree 



a = arg 

and penguin part. 

Take M.(B + — > K + K* ) in Eq. (1331) as an example, its "T" and "P" parts can be 
written as in the form of 



T 
P 



M aK Ct + F aK ( hjx + C 2 



(50) 



F e K 

+M eK 



1 1 1 

-Cs + C4 — -Cg — -C\o 
6 O 2 



^ r aK 



— ~Cq 

d 2 



C 5 - -C 7 



-C5 + Co + -C 7 + Cs 



+ M aK {Cs + C 9 



+M£ (C 5 + C r ) + FaK [ -Cs + C, + -Cg + C 10 



(51) 



In pQCD approach, the ratio z and the strong phase 5 can be calculated perturbatively. 
For B + — ► K + K*° and decays, for example, we find numerically that 



2.1, 



^(ir + K*°) 
2(iT* + K°) = 2.7, 



(5(K + K*°) 
5{K* + K°) 



-13°, 
-44°. 



(52) 



The major error of the ratio z and the strong phase 8 is induced by the uncertainty of 
Ub = 0.4 ± 0.04 GeV but is small in magnitude. The reason is that the errors induced by 
the uncertainties of input parameters are largely canceled in the ratio. 

From Eq. (JHJ), it is easy to write the decay amplitude for the corresponding charge 
conjugated decay mode 

M = V ub V* d T - V tb V t * d P = V ub V: d T [1 + ze«- a+ V] . (53) 

Therefore the CP-averaged branching ratio for B° — > KK* decay can be defined as 

Br = (\M\ 2 + |7W| 2 )/2 = \V ub V* d T\ 2 [l + 2^ cos « cos 5 + z 2 ] , (54) 

where the ratio z and the strong phase S have been defined in Eqs. (I48p and 
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CQ 



a(degree) 



FIG. 5: Branching ratios (in units of 10 
(dotted curve), B°/B° 
curve) as a function of CKM angle a. 



oiB+ -> K* + K° (dash-dotted curve), B+ 
(solid curve), B°/B° -> K+K*" + K~ K* + (dashed 



K K:* + K°K* '°^ A ^ ro'/— 



It is a little complicate for us to calculate the branch ratios of B°/B° — > /(/), since 
both 5° and 5° can decay into the final state / and / simultaneously. Due to B° — ~B° 
mixing, it is very difficult to distinguish B° from B . But it is easy to identify the final 
states. Therefore we sum up B°/B — > K°K* as one channel, and B°/B — > K K*° as 
another, although the summed up channels are not charge conjugate states 29]. Similarly, 

we have B°/B — > K + K*~ as one channel, and B°/B — ► K~K* + as another. We 
show the branching ratio of B°/B° -> K + K*~, B° /B° -> K~K* + , B + -> # + iT and 
5+ -> ir+lT decays clS cL function of a in Fig. [51 

Using the wave functions and the input parameters as specified previously, it is straight- 
forward to calculate the branching ratios for the four considered decays. The pQCD 
predictions for the branching ratios are the following 



Br(B 
Br(B 

Br(B°/B° 
Br(B°/B° -> K + K* 



<*0s 



K + K 
K* + K 



3.1±H(w b ) x 10 

h6.8 
-4.7 



-7 



18.3i5-?(w 6 ) x HT 7 , 



19.6l^(o; 6 ) x 10- 7 , 



5.4l 
1.0/ 



+ = 7.4^;° (u; 6 ) x 10~ 8 



(55) 
(56) 
(57) 
(58) 



where the major error is induced by the uncertainty of lo^ = 0.4 ± 0.04 GeV. 
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As a comparison, we also list the theoretical predictions in QCDF approach 



Br(B~ -»■ K~K*°) = 3.0±f;g x 10~ 7 , (59) 

Br(B~ -> iT~iO = 3.0i^;? x 10~ 7 , (60) 

Br(E° -> ZV°) = 2.6+^ x 10~ 7 , (61) 

Br(B° -> K°Z*°) = 2.9± 7 ;^ x 10~ 7 , (62) 

Br(B° -> K~K*+) = 1.4iJ°4 7 x 10" 8 , (63) 

Br(B° -> K + K*~) = 1.4+}° 4 7 x 10" 8 , (64) 



where the individual errors as given in Refs. |4( have been added in quadrature. For 
B~ — > K~K*° decay, the pQCD and QCDF predictions agree very well. For remaining 
decay modes, the pQCD predictions are larger than the QCDF predictions by a factor of 2 
to 5, although they are still consistent with each other within errors because the theoretical 
uncertainties are still very large. When compared with the experimental upper limits, the 
theoretical predictions in both approaches still agree with the data. The large differences 
between the pQCD and QCDF predictions will be tested by the forthcoming precision 
measurements. 



C. CP- violating asymmetries 

Now we turn to the evaluations of the CP-violating asymmetries of B — > KK* decays 
in the pQCD approach. For B + — ► K + K*° and B + — > K* + 7C° decays, the direct CP- 
violating asymmetries Afjp can be defined as: 

A dir = \M\ 2 -\M\ 2 = 2z sina sin 5 
CP |A4| 2 + |A4| 2 1 + 2zcosacos5 + z 2 ' 1 ' 

where the ratio z and the strong phase 5 have been defined in previous subsection and 
are calculable in PQCD approach. 

Using the definition in Eq. (I65j) . it is easy to calculate the direct CP- violating asym- 



metries for B ± —> ^K* (K*°) and B ± — > K* ± ~K"(K°) decays. The numerical results 



are 

1 dir I n± 



Aq P {B — > K K (K )) = -0.20 ± 0.05(a) ±0.02(w 6 ), 
A d j P (B ± -> K^K^K )) = -OA9t° °ol(a) ±0.07(cu b ). (66) 

for a = 100° ±20° and Uf, = 0.40 ±0.04 GeV. These pQCD predictions are also consistent 
with those in QCDF approach Q]: 

A d j P (B ± -> K+K^iK* )) = -0.24^;1, 

Aq^E^ — ► K* ±r K°(K )) = -0.13+°;^,. (67) 

where the individual errors as given in Ref. [4| have been added in quadrature. In Fig. [6l 
we show the a— dependence of the pQCD predictions of Afjp for — ► K ± K*° (K*°) (the 



solid curve) and B ± — > K*^K (K°) decay (the dotted curve), respectively. 
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FIG. 6: The direct CP asymmetry A^ p (in percentage) of B + — > K + K (the solid curve) and 
B+ -► K* + K (the dotted curve) as a function of CKM angle a. 

For B° I B° — » i^.fr^lK^.fr* ) decays, they do not exhibit CP violating asymmetry, 
since they involve only penguin contributions at the leading order, as can be seen from 
the decay amplitudes as given in Eqs. (1291 and ([3D]) . 

We now study the CP-violating asymmetries for B°/B — ► K + K*~(K~K* + ) decays. 
Since both 5° and 13° can decay to the final state K + K*~ and K* + K~ , there are four 
decay modes. Here we use the formulae as given in Ref. (29[. The four time-dependent 
decay widths for B°(t) -> K+K*-,B°(t) -> K-K*+,B°(t) -> and if(t) -> 

K+K*- can expressed by four basic matrix elements j29l]: 

<? = (K+K*-\H eff \B°), ft = (K + K*~ \H e ff\B°), 

g = (K-K* + \H eff \B^), h = (K-K* + \H eff \B°), (68) 

which determines the decay matrix elements of 5° — > K + K*~, B° — > K~K* + , B° — > 
K~K* + and I? — > K + K*~ at i = 0. The matrix elements g and ft, are given in Eqs. (1311) 
and (I32I) . The matrix elements ft and 7? are obtained from ft and g by simple replacements 
of £ u — > £* and £t — > i.e., changing the sign of the weak phases contained in the 
products of the CKM matrix elements £ u and £ t . 

Following the general procedure, the B° — B mixing can be defined as 

B 1 = p\B°) + g|S°), B 2 = p\B°) - (69) 



with |p| 2 + |g| 2 = 1. Following the notation of Ref. [29j, the four time-dependent decay 
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widths of the considered decay modes can be written as 

(\g\ 2 + \h\ 2 ) x {1 + <v cos(Amt) + a e+e , sin(Amt)} , 



T(B°(t) - 


■+ 


-K*~) 


= e" 


2 


T(B°(t) - 




~K*~) 


= e" 


2 


T(B°(t) - 


-> 1T 


-R* + ) 


= e" 


-r f l 
2 


T(B°(t) - 


-> if" 


-K*+) 


= e" 


-rtl 
2 



;(|<7| + W ) x {1 + cos(Amt) + a e+? > sin(Amt)} ,(70) 
where the four CP violating parameters are defined as 



M 2 -i/f . _-2M;f) 



g\* + \hr ^ l + \h/g\ 



2 



\W - \9\ 2 ~ 2 Mgf) .„. 

°-- = WTW'' ° ,+ "-TTW (71) 

where q/p = e 2l/3 . Using the decay amplitudes as given in Eqs. ( l31i) and ( 1321) . it is 
straightforward to calculate the above four CP-violation parameters. The central values 
of the pQCD predictions are 

a f > = 0.74, a e+t > = 0.68, 



a t = 0.25, a t+l > = -0.88, (72) 

for a = 100°. The a— dependence of these four CP violating parameters are shown in 
Fig. [3 It is difficult to measure these physical observables in current and forthcoming B 
meson experiments because of its tiny branching ratio ( ~ 10~ 8 ). 

At last, we will say a little more about the possible FSI effects. As mentioned in the 
introduction, we here do not consider the possible FSI effects on the branching ratios 
and CP-violating asymmetries of the B — > KK* decays. The FSI effect is in nature 
a subtle and complicated subject. The smallness of FSI effects has been put forward 
by Bjorken [30] based on the color transparency argument [&], and also supported by 
further renormalization group analysis of soft gluon exchanges among initial and final state 
mesons [2(j. At present, the excellent agreement between the pQCD predictions for the 
branching ratios and CP violating asymmetries and the precision measurements strongly 
support the assumption that the FSI effects for B — > Kir decays are not important 
[3]. For B — > KK decays, fortunately, good agreement between the pQCD predictions 



for the branching ratios of B + — ► K + K°, B° — ► K + K~ and K°K° decays jl5l . Il6| and 
currently available experimental measurements 3] indicates that the FSI effects are most 
possibly not important also [l6[]. Of course, more studies are needed about this issue, while 
further consistency check between the pQCD predictions and the precision data will reveal 
whether FSI effects are important or not. 



V. SUMMARY 

In this paper, we calculate the branching ratios and CP-violating asymmetries of 
-> K K*°(K°K* ), B°/B° -> K + K*-{K~K*+), B+ -> K + K*°, and B+ 
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FIG. 7: CP violating parameters of B°/B -> K + K*-(K-R* + ) decays: a e > (dash-dotted line), 
a? (dotted line) , a e+e i (dashed line) and a t+ i' (solid line) as a function of CKM angle a 

K* + K° decays, together with their charge-conjugated modes, by employing the pQCD 
factorization approach. 

From our calculations and phenomenological analysis, we found the following results: 

• The pQCD predictions for the form factors of B — > K and K* transitions are 

Ffr«{0) = 0.3518:S, A*- K * = 0A6+_° Z (73) 



for Ub = 0.40 ± 0.04 GeV, close to the light-cone QCD sum rule results [28 
the pQCD predictions for the CP-averaged branching ratios are 





Br(B + - 


> K + K*°) 


w 3.1 x 10 -7 , 




Br(B + - 


> K *+K°) 


w 18.3 x 10~ 7 


Br(B°/B° 


- K°K*° 


f K°K*°) 


w 19.6 x 10~ 7 


5r(5°/fi° - 


> K + K*- -f 


-K-R* + ) 


w 7.4 x 10" 8 . 



(74) 

The above pQCD predictions agree with the QCDF predictions within still large 
theoretical errors and close to currently available experimental upper limits. 

For the CP-violating asymmetries of the considered decay modes, the pQCD pre- 
dictions are generally large in magnitude. 
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APPENDIX A: NON-ZERO TRANSITION AMPLITUDES 

The factorizable amplitudes F eK *, F^„ and Ff^ t , F aK *, F^k* and F^ t have been 
given in Sec. IHII The remaining factorizable transition amplitudes in B — > KK* decays 
are written as: 

_ r°° 
F e x = A\f2nG pC p j 'K* m B / dx\dx 3 / bidbib 3 db 3 (pB{xi,bi) 

Jo Jo 

■ { [(1 + x 3 )<I>k(x3, h) +r K (l- 2x 3 ) {<p p K (x 3j b 3 ) + <f K (x 3 , 6 3 ))] 
■a s {t\)h e {x u x 3 , 6i, 6 3 ) exp[—S a (1^)] 

+2r K (p^(x 3 ,b 3 )a s (t 2 e )h e (x 3 ,xi,b 3 ,b 1 )exp[-S a (t 2 e )}} , (Al) 

*S = F eK, (A2) 

pi poo 

Ay/^TrGpCpfBrn^ / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

■ { [x 3 (j>K*(x 3 ,b 3 )(f)^ : (x2,b 2 ) + 2r K *r K <\) P K {x 2 ,b 2 ) ((1 + x 3 )(j) s K ,{x 3 ,b 3 ) 
-(1 - £ 3 )</4-,(a; 3 , 6 2 ))] a s (t 3 e )h a (x 2 , x 3 , b 2 , b 3 ) exp[-S d (t 3 e )} 

- [x2(f>K*{x3,b 3 )<f)K(x 2 ,b 2 ) +2r K »r K (/) s K ,(x 3 ,b 3 ) ((1 + x 2 )<Pk{ x 2, b 2 ) 

-(1 - x 2 )(f)K(x 2 ,b 2 ))] a„(tf)h a (x 3 , x 2 ,b 3 ,b 2 ) exp[-S d (t*)]} , (A3) 

FaK = 'FaK, (A4) 

rl POO 

F P K = %\f2G f kC pm B f b / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

■ { [2r K (j) K * (x 3 , 6 3 )0^(x 2 , b 2 ) + x 3 r K * (<p s K * (#3, 63) - 0k* O3, ^2)) </4 O2, 62)] 
•a s (tg)/i a (x 2 , x 3 , b 2 , b 3 ) exp[-S d (tl)} 

+ [2r K *<f) s K * (x 3 , b 3 )(p^(x 2 , b 2 ) + x 2 r K (<f> K {x 2 , b 2 ) - 0^(x 2 , b 2 )) (j) K * (x 3 , b 3 )] 
■a s (ti)h a (x 3 ,x 2 ,b 3 ,b 2 )exp[-S d (tj)]} . (A5) 

For B — > KK* decays, the non-factorizable transition amplitudes not shown explicitly 
in Sec. IHII are written as: 

16 

M e K = -=GpTrCpm B / dx\dx 2 dx 3 I bidbib 2 db 2 (j)B(xi,bi)(f)K*(x 2 ,b 2 ) 

v3 Jo Jo 

■ { [-x 2 <Pk(x 3 , b 2 ) + r K x 3 (<j> K (x 3 , b 2 ) - (j) K (x 3 , b 2 ))] 

■a s (t))h){x x , x 2 , x 3 , bi, b 2 ) expl-S^t})] 

- [(x 2 -x 3 - l)(j)^(x 3 , b 2 ) + r K x 3 ((f) K (x 3 , b 2 ) + (j) K {x 3 , b 2 ))] 

■a s (tj)hj(x u x 2 ,x 3 ,b 1 ,b 2 )exp[-S b (t 2 f ))} , (A6) 
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16 f 1 [°° 

= — -=G F irC F r K *m B / dx x dx 2 dx 3 I b 1 db 1 b 2 db 2 (f) B (x 1 ,bi) 
v3 Jo Jo 

■ { [x 2 <t>K{x-i, h) ((f) s K »(x 2 ,b 2 ) - 0^(x2,6 2 )) -r K (x x (<f)^(x 3 ,b 2 ) 
-(t> T K {x 3 ,b 2 )) {<t>K*{x2,b 2 ) - 0^.(2:2,62)) -^2 (0^(^3,62) 
-4>K(x 3 ,b 2 )) (4> s K *(x 2 , b 2 ) - <j> t K .{x2,b 1 )) ~ x 3 ((j)K(x3,b 2 ) + <pK( x 3,b 2 )) 

■ ((f) s K *{x 2 ,b 2 ) +(f> t K *(x 2 ,b 2 )))] a a (i})/i}(a;i,X2,X3,6i,62)exp[-5' 6 (i})] 
+ [(1 - x 2 )0^(x 3 , b 2 ) [<j> s K , (x 2 , b 2 ) - <\) l K , (x 2 , b 2 )) - r K [x\ (0£(z 3 , h) 

-4>K( x 3,b 2 )) (<j) s K *( x 2,b 2 ) -^.(x 2 ,6 2 )) - (1-^2) (0^(2:3,62) 
-0^(x 3 , 62)) (ar 2 , 62) - <f>K* ( x 2, b 2 )) - x 3 ((/>k(x 3 , b 2 ) + 0^(x 3 , b 2 )) 

■ {<t>K*{ x 2, b 2 ) + <j> t K .(x 2 ,b 2 )))] a s (^)/iJ(x 1 ,x 2 ,X3,6 1 ,6 2 )exp[- ( 5 b (^)]} , (A7) 

16 Z" 1 /"°° 

M aK = -=nG F C F m B / dx ± dx 2 dx 3 I b 1 dbib 2 db 2 4>b{x\, h) 

v3 Jo Jo 

■ { [ x 2<Pk* ( x 2, h)(p^(x 3 , b 2 ) + r K *r K (((x 2 + x 3 + 2) 0^, (x 2 , 6 2 ) 

+ (x 2 - x 3 ) 0^» (x 2 , b 2 )) 0^(x 3 , 6 2 ) + 0^(a; 3 , b 2 ) (-x 3 (0^* (x 2 , 62) 
-0^(x 2 ,6 2 )) -20' x *(x2,6 2 ) + x 2 (<f) s K *(x 2 ,b 2 ) + 0' x *(x 2 , 62))))] 
• a s (tj)/ij(xi, x 2 , x 3 , 61, b 2 ) exp[-S c (t 3 f )] 

- [ x 3<Pk* ( x 2, b 2 )(p^(x 3 , b 2 ) + r K «r K [x 2 (0^(x 3 , b 2 ) - <Pk(x 3 , b 2 )) 

■ (<p s K *( x 2,b 2 ) -4>k* ( x 2,h)) +x 3 (0£(x 3 , b 2 ) + 4>k(x 3 , b 2 )) 

■ ((f) s Kr (x 2 ,b 2 ) + (/> t K .(x 2 ,b 2 )))] a s (tj)hj(x u x 2 ,x 3 ,b 1: b 2 )ex.p[-S c (tj)]} ,(A8) 
Kk = M*., (A9) 
<1 = (MO) 

where = /m B with = m 2 K /{m s + m d ). 



APPENDIX B: RELATED FUNCTIONS 

We show here the function hiS, coming from the Fourier transformations of H^°\ 

h e (x 1 ,x 3 , 61, 63) = Ko {y/x^rriBbi) [0(&i - 63)^0 (v^3 m iA) ^0 {V^^sh) 

+6(b 3 - b^Ko (^m B b 3 ) Iq (V^mBbi)} S t (x 3 ), (Bl) 

^0(^2, x 3, h, h) = K (iy/x 2 x 3 m B b 2 ) [9(b 3 - b 2 )K (i^fxl,m B b 3 ) I {iyfxl,m B b 2 ) 

+6(b 2 - b 3 )K {i^x~ 3 m B b 2 ) Io {i^ 3 m B b 3 )\ S t (x 3 ), (B2) 

,x 2 ,x 3 ,b 1 ,b 2 ) = j#(& 2 - bi)l (M By /xix 3 b 1 )K (M Bv /xix 3 b 2 ) 

} f K (M B D (j) b 2 ), for D 2 {j) >0\ 
+ ~ h) \ ■ [j^\M B ^\ b 2 ), for Df, < oj ' ™ 

22 



h f (x 1 , x 2 , x 3 , 61, b 2 ) = < 6{pi - b 2 )K (i^X2X 3 b 1 M B )l (i^X2X3b 2 MB) + (h <-+ b 2 ) 



•—U^\s/xx + x 2 + x 3 - xix 3 - x 2 x 3 bi M B ) 



(B4) 



h 4 f (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 



+ (b 



9(bx - b^Koii^x^sbiMB^iy/x^foMB 
K {M B F {1) bx), 



for F { 2 1} > N 



^\M B J\F^\b 1 ), for F ( 2 1)<0 



, (B5) 



where j=l and 2, Jo is the Bessel function and Kq, Iq are modified Bessel functions 
K (—ix) = —(ir/2)Y (x) + i(ir/2)J (x), and F?^ , D/j^s are defined by 



n 2 
n 2 



(i) ' ^U) 

{x x - x 2 )x 3 , 
(xx - x 2 )x 3 , 
-(1 - Xx - x 2 )x 3 . 



The threshold resummation form factor S t (xi) is adopted from Ref . [24 

2 1+2 T(3/2 + c 



S t (x) 



VtFT(1 + c) 



[ar(l -x)] c 



where the parameter c = 0.3. This function is normalized to unity. 
The Sudakov factors used in the text are defined as 

S a (t) = s (xxm B / V2, bx^j + s (x 3 m B / y/2, 63) +s((l- x 3 )m B / y/2, b 3 J 



1 

Ti 



In 



ln(*/A) 



In 



ln(t/A) 



S b (t) 



ln(6iA) -ln(M-). 
s {xiItlb/ V2, bxj + s (x 2 m B / yp2, b 2 ^j + s ^(1 - x 2 )m B / y/2, b 2 ^j 

+s (x 3 m B / y/2, bxj + s ((1 - x 3 )m B / y/2, h 
1 



ln(6iA) ' *"-ln(6 2 A)_ 
s \ Xxm B / y/2, bxj + s (x 2 m B /y/2, b 2 ) + s ( (1 - x 2 )m B / y/2, b 2 ) 

+s (x 3 m B / y/2, 6 2 ) + s f (1 - x 3 )m B /V2, b 2 



1 



In 



ln(t/A) 



2 In 



ln(t/A) 



S d (t) 



\n(bxA) -ln(62 A) J ' 

= s (x 2 rn B / y/2, b 2 J + s (x 3 m B / y/2, b 3 J +s(jl — x 2 )m B / y/2, b 2 J 

+ s ((1 - x 3 )m B /V2, 6 3 ) - 



In ln(tM) , I In ln(tM) 



ln(6 2 A) 



ln(6 3 A) 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 



(Bll) 
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where the function s(q, b) are defined in the Appendix A of Ref. 22]. The scale Ws in the 
above equations are chosen as 

t\ = max( v / a^m B , 1/6 X , l/6 3 ) , 
t\ = max( v /a^"m jB , 1/63) , 



tr e = max( v /i^m B , l/6 2 , l/k 



■3 



■3 J 



t e = max( v /x 2 ~m B , l/6 2 , 1/fe; 

= max( v /i^m B , a/O^i _ x 2 )x 3 m B , l/&i, l/6 2 ) , 

^ = max( v /xiX 3 m B , a/(1 - Xi - x 2 )x 3 m B , I/61, l/6 2 ) , 

= max( v / a;i + x 2 + x 3 - - x 2 x 3 m B , \Jx 2 x 3 m B , l/&i, 1/&2) , 

= max( v / (xi - x 2 )x 3 m B , y/x 2 x 3 m B , l/b u l/b 2 ) . (B12) 
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